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Primordial Black Holes (PBHs) are black holes that may have been created in the early Universe 
and could be as large as supermassive black holes or as small as the Planck scale. It is believed that 
a black hole has a temperature inversely proportional to its mass and will thermally emit all species 
of fundamental particles. PBHs with initial masses of ~ 5.0 x 10^“* g should be expiring today 
with bursts of high-energy gamma radiation in the GeV/TeV energy range. The High Altitude 
Water Cherenkov (HAWC) observatory is sensitive to the high end of the PBH gamma-ray burst 
spectrum. Due to its large field of view, duty cycle above 90% and sensitivity up to 100 TeV, the 
HAWC observatory is well suited to perform a search for PBH bursts. We report that if the PBH 
explodes within 0.25 light years from Earth and within 26 degrees of zenith, HAWC will have a 
95% probability of detecting the PBH burst at the 5 sigma level. Conversely, a null detection from 
a 2 year or longer HAWC search will set PBH upper limits which are significantly better than the 
upper limits set by any previous PBH search. 
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1. Introduction 

In the current Universe, there are no known processes that can create black holes with masses 
less than a few solar masses. However, early in the Universe, conditions may have been extreme 
enough to create black holes with masses ranging from the Planck mass to super-massive black 
holes [1]. These black holes are called Primordial Black Holes (PBHs). PBH production in the 
early Universe can have observable consequences, spanning from the largest scales, for example 
influencing the development of large-scale structure in the Universe, to the smallest scales, for 
example enhancing local dark matter clustering. In the present Universe, PBHs in certain mass 
ranges may constitute a fraction of the dark matter [1]. 

It is expected that every black hole thermally radiates (‘evaporates’) with a temperature in¬ 
versely proportional to its mass [2]. The emitted radiation consists of all fundamental particle 
species with masses less than about the black hole temperature [3]. PBHs with an initial mass of 
~ 5.0 X 10*“^ g should be expiring now with bursts of high-energy particles, producing gamma ra¬ 
diation in the GeV - TeV energy range [4]. Positive detection of a PBH burst event would provide 
important information about the processes that shaped the early Universe. In addition, confirmed 
detection will give access to particle physics at energies higher than currently achievable by ter¬ 
restrial accelerators. Even the non-detection of PBH burst events in dedicated searches would give 
important information about fundamental physics and the early Universe. 

Numerous observatories have searched for PBH burst events using direct and indirect methods. 
These methods are sensitive to the PBH distribution at various distance scales. Observatories that 
observe photons or antiprotons at ~ 100 MeV can probe the cosmologically-averaged or Galactic- 
averaged PBH distribution whereas TeV observatories directly probe PBH bursts on parsec scales. 
Due to the possibility that PBHs may be clustered at various scales, all these searches provide 
important information. We also note that the TeV direct search limits apply to any local bursting 
black holes, regardless of their formation mechanism or formation epoch. Table 1 gives a summary 
of various search methods, the distance scales they probe and their current best limits. 


Distance Scale 

Eimit 

Method 

Cosmological Scale 

< lO^^pc^^yr^^ 

(1) 

Galactic Scale 

< 0.42 pc^^yr^^ 

(2) 

Kiloparsec Scale 

< 1.2 X 10^^ pc^^yr^' 

( 3 ) 

Parsec Scale 

< 1.4 X 10“^ pc^^yr^' 

( 4 ) 


Table 1: PBH burst limits on various distance scales; (1) from the 100 MeV extragalactic gamma-ray 
background assuming no PBH clustering [1, 5], (2) from the Galactic 100 MeV anisotropy measurement [6], 
(3) from the Galactic antiproton flux [7] and (4) from Very High Energy direct burst searches [8]. 


In this paper, we present the sensitivity of the High Altitude Water Cherenkov (HAWC) obser¬ 
vatory to PBH bursts, which tries to improve the PBH limit at the parsec scale. 

2. HAWC Observatory 

The HAWC observatory is a very high energy (VHE) observatory located at Sierra Negra, 
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Mexico at an altitude of 4,100 m above sea level. HAWC consists of 300 water tanks, each 7.3 m 
in diameter and 4.5 m deep. Each tank houses three 8-inch photo multiplier tubes (PMTs) and one 
10-inch PMT. The PMTs detect Cherenkov light from secondary particles created in extensive air 
showers induced by VHE gamma rays of energies between ~50 GeV and 100 TeV. The main data 
acquisition system measures the arrival direction and energy of the VHE gamma ray by timing the 
arrival of the secondary particles on the ground and measuring the amplitude of the PMT signals. 
The direction of the primary particle may be resolvable with an error of between 0.1 and 2.0 degrees 
depending on its energy and location in the sky. HAWC has a large instantaneous field-of-view of ~ 
2 sr and will have a high duty cycle of greater than 95%. Thus, HAWC will be able to observe high- 
energy emission from gamma-ray transients. Eor more information about the HAWC observatory 
please see references [9, 10] and references therein. 

3. Methodology 

3.1 Primordial Black Hole Burst Spectrum 

As the PBH radiates, it continually loses mass and its temperature increases to very high en¬ 
ergies. The manner in which the PBH expires depends on physics at very high energies. In the 
Standard Evaporation Model (SEM) [3, 11] (based on the Standard Model of particle physics), a 
PBH should directly radiate those fundamental particles whose Compton wavelengths are com¬ 
parable with the size of the black hole. When the black hole temperature exceeds the Quantum 
Chromodynamics (QCD) confinement scale (~250-300 MeV), quarks and gluons should be radi¬ 
ated [3, 5]. The quarks and gluons will then fragment and hadronize as they stream away from the 
black hole, analogous to the jets seen in terrestrial accelerators [3, 4]. On astrophysical timescales, 
the final jef producfs will be photons, neufrinos, elecfrons, posifrons, profons and anfi-profons. In 
fhis work, we assume fhe SEM as our emission and particle physics model. Please refer fo refer¬ 
ence [12] for more information abouf fhis model, and for fhe calculafion of PBH bursf lighf curves 
in fhe HAWC energy range and fheir energy dependence. 

In fhe SEM model, fhe black hole femperafure (T) can be expressed in ferms of fhe remaining 
evaporafion lifefime (t) of fhe black hole (fhaf is, fhe lime lefl until fhe black hole stops evaporating) 
as follows [11, 13]: 



(3.1) 


for lemperalures well below ils ullimale or Planck femperafure. The emission rale increases as fhe 
black hole shrinks [5]. Eor black holes wilh lemperalures greafer lhan several GeV al fhe slarl of 
fhe observalion, fhe lime inlegraled pholon flux can be paramelerized as [13] 



(3.2) 


for gamma-ray energies E > 10 GeV. This paramelrizalion includes bolh fhe direclly radialed pho¬ 
tons and Ihose produced by fhe decay of olher direclly radialed species. Eigure 1 shows fhe lolal 


3 









Sensitivity ofHAWC to PBH Bursts 


J. H. MacGibbon 



Figure 1: Time-integrated gamma-ray spectrum over various PBH remaining lifetimes using the 
parametrization of Equation 3.2. The black hole temperature at the start of observation is shown in paren¬ 
theses. 


gamma-ray spectrum for various PBH remaining lifetimes ranging from 0.001 s to 100 s. The 
parametrization of the instantaneous PBH spectra is presented in reference [12]. 


3.2 PBH Upper Limit Estimation 

In the case of a null detection by HAWC, we can derive an upper limit on the local PBH 
burst rate density. The parametrization given in equation 3.2 can be used to calculate the expected 
number of photons detectable by an observatory on the Earth’s surface. For a PBH burst of duration 
T at a non-cosmological distance r and zenith angle 6, the number of expected photons is 


n{r,e,z) 


( 1 -/) 

4Kr'^ 


/■^2 dN{T) 
'e, dE 


A{E,e)dE 


(3.3) 


where dN/dE is the PBH gamma-ray spectrum integrated from remaining lifetime T to 0. The 
energies Ey and E 2 correspond to the energy range of the detector, A{E,d) is the effective area of 
the detector as a function of photon energy and zenith angle, and / is the dead time fraction of the 
detector. 

The minimum number of counts needed for a detection for different burst durations T, /To (0;, t) , 
can be estimated by finding the number of counts required over the background for a 5a signifi¬ 
cance. We define a 5a defecfion affer correction for Nt trials to be the number of counts n which 
would have a Poisson probability P corresponding to a p-value pc given by 


Pc = Po/Nt = P{> n\nbk) 


(3.4) 
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where po = 2.3 x 10^’ (the p-value corresponding to 5a), nbk is the number of background counts 
expected over the burst duration T, and P{> n|nbk) denotes the Poisson probability of getting n or 
more counts when the Poisson mean is nbk- We take the value of ;Uo(0,',t) to be the amount of 
expected signal which would satisfy this criterion 50% of the time. To find iJ.o{6i,x), we estimate 
the value at which the Poisson probability of finding at least n counts is 50% by using the relation 

/’(>«|(«bk + Ato(0;,T))) =0.5. (3.5) 


The maximum distance from which a PBH burst could be detected by HAWC can be calculated by 
equating the values of /To (0,, t) to p {r, 0,-, t) and solving for r, 


rmM,r) = .p , T' ^^A{E,ei)dE. 


4Kpo{di,T:) Jei dE 


(3.6) 


The detectable volume from which a PBH could be detected is then 


KW = Ev(e,.i) = 


(3.7) 


where FOVeff(6i) denotes the effective field-of-view of the detector for a given zenith band i, 


FOVeff(0;) — 27r(cOS 0,’ fnin COS0; fnax) Sr 


(3.8) 


and 0,,rnin and 0,,max are the minimum and maximum zenith angles in band i. 

If we assume that the PBHs are uniformly distributed in the solar neighborhood, the X% 
confidence level upper limif (ULx) on fhe local rale densily of PBHs bursls (fhal is, Ihe number of 
bursls occurring locally per unil volume per unil lime) can be eslimaled as 

m 

ULx = —, (3.9) 

if, al fhe X% confidence level, fhe defector observes zero bursls. Here V is fhe effeclive deleclable 
volume from which a PBH can be delecled, S is fhe lolal search duration and m is fhe upper limif on 
fhe expecled number of PBH bursls given fhal al fhe X% confidence level zero bursls are observed 

al Earth. For Poisson flucluations, m = ln(l/(l —X)). If X = 99%, Ihe upper limit on the PBH 

burst rate density will be 

4.6 

ULg9 = —. (3.10) 


4. Results and Discussion 

The methodology presented in Section 3 can be used to find Ihe projected PBH bursl rate 
densily upper limils for HAWC in Ihe case of null detection in a 5 year search. We employ a Monte 
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Figure 2; PBH Burst Rate Density Upper Limits projected for HAWC [14], compared with limits from 
previous direct search experiments. 

Carlo simulation which models the interaction of photons and cosmic rays with the atmosphere and 
the response of the detector to the extensive air showers they generate. The effective area A{E,6) 
is equated to the ratio of the number of events that satisfies a given set of cuts to the total number of 
events multiplied by the total throw area of the Monte Carlo simulation. The cuts are comprised of 
a trigger cut, an angle cut and a gamma-hadron separation cut. For the trigger cut, HAWC accepts 
events with the number of PMTs hit by the air shower, nHit, greater than a certain value. The angle 
cut is used to specify the direction of the photons and is a measure of HAWC’s angular resolution. 
In the simulated events, an angular parameter defined fo be fhe difference befween fhe frue location 
of fhe particle in fhe sky and fhe reconsfrucfed sky locafion of fhe parficle, DelAngle, is used as a 
proxy for fhe angular search bin-size. 

Because we seek fhe sensifivify in fhe case where fhere is no prior knowledge of fhe bursf 
locafion, we need fo fake info accounf fhe number of frials performed for fhe search. The optimum 
spafial bin-size depends on fhe search durafion, fhe trigger criteria, and the value of the gamma- 
hadron separation parameter. The number of time bins is estimated by dividing the total search 
period (estimated as 5 years for HAWC) by the burst duration T. Thus the total number of trials 
depends on the burst duration T, the optimal spatial bin-size DelAngle, the trigger criterion nHit 
and the value of the compactness parameter. In order to find fhe optimum sef of cufs, we performed 
a simple parameter search and idenfified fhe sef of values which give fhe sfrongesf PBH bursf rate 
upper limif. We refer fhe reader fo reference [14] for more defails on fhe opfimizafion process and 
fhe background rale calculations. 

Figure 2 displays fhe PBH Bursf Rale Density Upper Limifs projecled for HAWC, compared 
wilh limifs from previous direcl search experimenfs. The HAWC projecled PBH rale density limit 
is strictest around 10 seconds. The general features of Figure 2, which shows limits as a function 
of T, can be understood as follows. From Equations 3.6-3.10, we can see that UL, the upper limit 
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on the local PBH rate density, scales as 

ULr. (4.1) 

where /Tq is the sensitivity for a given search bin T, and Ny{T) is the number of observable photons 
produced by the source. Better sensitivity corresponds to smaller /Tq, the number of signal photons 
required for detection of a signal, and a stricter PBH limit. For a source at a given distance (that is, 
at the outer edge of the volume considered), Nj{x) decreases when the search interval is smaller, 
and produces a weaker PBH limit. Shorter intervals incur less background events but fewer source 
photons. In this case, /Xo is dominated by statistical fluctuations, in particular those associated 
with the detector background rate, while Ny{'t:) is dominated by the PBH emission time profile, 
slightly modified by the energy-dependent effective area of the detector. These dependencies are 
both power laws, as seen in Figure 1 of reference [12], and, despite their very different physical 
origins, nearly cancel. Secondary effects such as the larger number of trials incurred for shorter 
intervals (proportional to 1 /t), and the ability to optimize background rates for larger T for which 
the detection efficiency is higher, give the residual T dependence seen in our Figure 2. 

According to our study, if a single PBH explodes within 0.074 parsecs of Earth and within 26 
degrees of zenith, HAWC will have a 95% probability of detecting it at 5 a (as optimized in a 10 s 
search after trials corrections). Also, HAWC would see with 95% probability a PBH burst within 
37 degrees of zenith if it happens within 0.058 parsec of Earth [14]. 

The HAWC estimated sensitivity presented above was based on choosing a single set of se¬ 
lection criteria (and gamma-hadron selection) for each search window duration. We have begun 
exploring improvements which might be obtained by making these selection criteria depend on 
nHit, the number of PMTs participating in the event, as is used in current HAWC analyses [9]. We 
evaluated HAWC’s PBH sensitivity using -dependent cuts tuned for these other HAWC analy¬ 
ses and see a factor of 1.3 improvement in the expected PBH burst-rate-density limit for 1 and 10 s 
search durations. Once the cut optimization is performed for a full likelihood search incorporating 
knowledge of the PBH time and energy dependence discussed in reference [12], we may obtain a 
larger improvement. 
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